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Abstract 

We construct and discuss solutions of 50(1, 2) x 5*0(1, 2) Chern-Simons theory which correspond to 
multiple BTZ black holes. These solutions typically have additional singularities, the simplest cases 
being special conical singularities with a 2ir surplus angle. There are solutions with singularities 
inside a common outer horizon, and other solutions with naked conical singularities. Previously 
such singularities have been ruled out on physical grounds, because they do not obey the geodesic 
equation. We find however that the Chern-Simons gauge symmetry may be used to locate all 
such singularities to the horizons, where they necessarily follow geodesies. We are therefore led to 
conclude that these singular solutions correspond to physically sensible geometries. 

Boundary charges at infinity are only sensitive to the total mass and spin of the black holes, and 
not to the distribution among the black holes. We therefore argue that a holographic description in 
terms of a boundary conformal field theory should represent both single and multiple BTZ solutions 
with the same asymptotic charges. Then sectors with multiple black holes would contribute to the 
black hole entropy calculated from a boundary CFT. 
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1 Introduction 



Three-dimensional gravity has been a useful laboratory for exploring quantum gravity in a sim- 
plified setting. For a negative cosmological constant there are black hole solutions 0, and the 
Bekenstein-Hawking entropy of these BTZ black holes have been attributed to boundary degrees 
of freedom at the horizon J3|, U , at infinity |J or at any intermediate timelike surfaces || . 

Strominger's asymptotic approach makes use of a particular property of asymptotically AdS 
solutions of three-dimensional gravity discovered by Brown and Henneaux jffl: the asymptotic 
isometries are represented canonically by a Virasoro algebra. The BTZ mass and spin determine the 
transformation properties under the conformal transformations (conformal weights) . In conformal 
field theory an argument by Cardy Q can be used to relate central charges of Virasoro algebras 
to the densities of states at high weights. Similarly the asymptotic density of states of quantum 
gravity is fixed by the central charge. It is found to agree with expectations from the BTZ horizon 
area. This elegant argument is independent of the precise conformal field theory representing 
quantum gravity, and gives few details about the theory. Carlip has combined it with his horizon 
approach, to open the way to an understanding of the universal nature of black hole entropy |^|, 
and its relation to horizon area. The price is that the horizon is treated as an input, rather than 
as a consequence of the global geometry, and again that details of the field theory are lost. 

The ability to compute black hole entropy quantum mechanically does not mean that it is 
fully explained. Even in the simple case of 2+1 dimensions the gravitational backgrounds that 
contribute to the entropy are poorly understood. Ideally a correct count of boundary degrees of 
freedom at the horizon or at infinity should also tell us what bulk geometries are relevant, and how 
they are excited. They may also be represented differently in different quantum gravity theories. 
(In 2+1 dimensions there are several inequivalent quantizations pp|.) To start investigating what 
geometries may represent the entropy we choose to study one description, Chern-Simons theory 
p+| , fL2[ . In Chcrn-Simon theory the map to the boundary theory is well-known and produces a 
WZW theory gg,[l|]. 

On the classical level one could ask what constant curvature metrics (solving the equations of 
motion) look asymptotically like BTZ black holes, and could be expected to be equally important 
as the standard BTZ solution for the black hole entropy. Bahados Jl6| (see also [[TtJ and jisj]) 
has given a simple analytic and general characterization of such solutions, but unfortunately the 
analytic expression of the solution does not give directly the geometric structure of the spacetime. 
In string theory approaches to black hole entropy, BPS solutions which can be separated into multi- 
source solutions play a prominent role [fl9f . This indicates that similar solutions may be of interest 
also in pure gravity. Indeed, we will find that asymptotically, Chern-Simons multi-source solutions 
typically are closer than most of the solutions in |l^] and jl7) to the standard BTZ solutions. 

Chern-Simons multi-source solutions have been discussed as candidates for stationary multi- 
black hole solutions by Coussaert and Henneaux pOfl . Clement [^lj found similar solutions in 
a metrical formulation and generalized them to dynamical solutions with moving sources. His 
main motivation for doing so was to correct the shortcoming also observed in pQ] , that the sta- 
tionary solutions necessarily involve additional conical singularities which typically do not follow 
geodesies^]. A non-geodesic behaviour signals an unwanted transport of energy-momentum between 
the singularity and spacetime. 

From our perspective this state of affairs is quite puzzling. The multi-source solutions could be 
expected to be on equal footing with the BTZ solutions on the basis of their asymptotic behaviour, 
and as Chern-Simons solutions they are no less regular. Hence calculations of black hole entropy 
based on Chern-Simons theory and its expression in terms of a boundary WZW model, which 
reproduce Bekenstein-Hawking's result, appear to include unphysical geometries. Is Chern-Simons 
theory, which has attracted so much attention as a model of quantum gravity, just not sensible as 

lr The presence of these singularities indicates that such solutions are different from the multi-black hole solutions 
with mul t i ple asymptotic regions that have been discussed by Brill [E3, El and the wormhole solutions by Bengtsson 
et al. |i[ P 
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a theory of gravity? 

Fortunately there is a caveat in the above argument. We find that multi-source solutions have 
similar asymptotic behaviour, but also that they represent physical (though singular) geometries. 
The crucial observation is that one should deal with gauge equivalence classes of solutions to Chern- 



Simons theory rather than with individual solutions. As we show in section 4.2.1, all "unacceptable" 
Coussaert-Henneaux solutions are gauge equivalent to perfectly acceptable solutions^. This is 
possible because the Chern-Simons gauge group is larger than the diffeomorphism group. For a 
discussion of how the gauge symmetry conspires with the presence of degenerate metrics in the 
Chern-Simons 'formulation' see Matschull p^j. In practice, gauge transformations move the static 
conical singularities to the horizons, where they obey the geodesic equation (by infinite redshift) 
as already observed by Clement. If all Chern-Simons solutions have similar sensible representative 
geometries or not is left as an open question. The Coussaert-Henneaux solutions dealt with here 
certainly constitute an important subclass. 

In section ^ we give our Chern-Simons formulation of the BTZ black hole and in section || we 
write this Chern-Simons solution in a more geometric way, and are led to a much more general 
solution, which includes the multi-black hole solutions, some of which have already been discussed 
by Coussaert and Henneaux. We also discuss how gauge transformations act on all these these 
solutions. In section [| we specialize to the case of two black holes (more precisely two excluded 
regions with closed timelike curves). We study the properties of this solution and the role of 
degenerate metrics. We solve the problem of the non-geodesic singularities of the Coussaert- 



Henneaux solutions by choosing a suitable gauge in 4.2.1 and we end with conclusions in section 
I 

2 The BTZ black hole 

In 2 + 1 dimensional gravity with a negative cosmological constant there exists a black hole solution 
to Einstein's equations, the BTZ black hole, jlj]. It can be viewed either as a metric approaching 
an AdS form asymptotically, or as a quotient of Anti-de Sitter space Q . The BTZ- metric can be 
written, 

ds 2 = ~N 2 dt 2 + N- 2 dr 2 + r 2 {N^dt + defy 2 (1) 
where the lapse function N and the angular shift are 

N<=' J=^ (3) 

and < r < oo, — oo < t < oo, < 4> < 2ir. M is the mass of the black hole and J is the angular 
momentum. Both these quantities can be expressed in terms of the values of r (r + and r_ ) when 
the lapse function N vanishes. They correspond to the outer (r + ) and the inner (r_) horizon of 
the black hole. For the horizons to exist we need M > and | J| < Ml. When r + coincides with 
r_, we get extremal black holes, |J| = Ml. We will be concerned mainly with the non-extremal 
case. 

The cosmological constant A is related to the length scale I by the A = —111 2 . We choose 



units such that 1=1. To facilitate the Chern-Simons formulation in section 2.2 wc rewrite the 



metric differently the outer region r>r + , the intermediate region r + > r > r_ and the inner region 

2 It may seem strange that some gauge potentials give meaningful metrics while others in the same class do 
not, but this property is in fact intrinsic to the Chern-Simons approach. The solutions to the equations of motion 
are pure gauge, and they may locally be transformed away, giving a completely degenerate metric, unless further 
conditions on the vector potentials are imposed. 
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r_ > r > 0. Thus we make the following Rindler-likc coordinate transformation in each region 
r > r+, r + > r > r_ and r_ > r > 0: 



I : r 2 — r\ cosh 2 (p — a — — ) — r 2 sinh 2 (p — a — — ), a + — < p < oo (4) 



II: r 2 = r 2 _ cos 2 (p - a) + r\ sin 2 (p - a), a<p<a+^ (5) 
III: r 2 = r 2 _ cosh 2 (p - a) - r\ sinh 2 (p - a), 0<p<a (6) 



a — arctanh ( - — J (7) 

The constant a is choosen in such a way that r = corresponds to p = 0. In these coordinates we 
get a one to one correspondence between r and p. This will lead to the following metrics: 

I : ds 2 — — sinh 2 ( ( o — a — — ) [r+dt — r-d(ft\ 2 + dp 2 
+ cosh 2 (p — a — — ) [r_G?t — r + d<ft} 2 

II : ds 2 = sm 2 (p - a) [r_di - r + d0] 2 - dp 2 (g) 
+ cos 2 (/9 - a) [r+eft - r-d<ft\ 2 
III : ds 2 = - sinh 2 ( j o - a) [r-dt - r+dcf)] 2 + dp 2 
+ cosh 2 (p — a) [r+dt — r-dcf)] 2 

If we look at the metric in the inner region III we find that our choice of a causes the coefficient 
of d(j) 2 to vanish precisely when p = and to become negative when p < 0, i.e. we will get closed 
timelike curves (CTCs). Excluding the negative p region in fact removes all CTCs ||. We also 
note that t is always a global Killing coordinate, timelike in / and spacelike in II and III. The 
xy plane is euclidean in J, lorentzian in II and euclidean in iTJ, implying that light cones are 
drastically tilted inside the black hole. The radial coordinate p is spacelikc in /, timelike in 77 and 
spacelike in 777. 

Here it makes sense to pause and think about the split into three different coordinate regions. 
The point we want to make may seem trivial in the metric formulation, but it will reappear in the 
Chern-Simons formulation. Although the boundaries between the regions happen to coincide with 
the positions of the inner and outer horizons there is of course nothing special going on locally in 
these places. So why do we not simply continue our expressions from one side of the boundary to 
the other instead of changing analytic forms from region to region? The answer is that the analytic 
expressions of the metric (^) become degenerate at the boundaries of the regions, indicating that 
the coordinates become singular there. In fact, if we were to use the region I expression for all p 
we would still have a spacetime divided into two separate regions because Einstein's equations of 
motion cannot really be applied to this degenerate metric. The true rationale for the matching 
of different metrics across the boundaries between regions is that one can find a coordinate chart 
covering the boundary, and diffeomorphic transformations on either side to the respective forms of 
the metric. 



2.1 Chern-Simons formulation of gravity 

In the Chern-Simons formulation of three-dimensional gravity [ fill isometries of the AdS back- 
ground are gauged. For AdS the isometry group is SO(l, 2) x 5*0(1, 2) and we call the respective 
gauge fields of each factor A = A k Jk and A — A k 3k- The SO(l, 2) generators of a factor 
of the group are different from those of the other factor, but since they never appear multiplied 
together we shall not distinguish between them. The commutation rules within each factor are 
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[Jfe, Ji] = e m kiJ m , with the convention e°i 2 = — £012 = — 1, and metric r\ a h = 2Tr(J a Jb) of signature 
(—1,1,1). The Chern-Simons three- form 

Trj^ AdA + ^AAA Aa| (9) 

and its counterpart for the other factor then serve as Lagrangian densities, which automatically 
yield a generally covariant action. The equations of motion 

F = dA + AAA = and F = dA + AAA = (10) 

are then actually equivalent with Einstein's equations, provided the identifications 

9ij = e?e b 3 Vab (11) 

e? = \(A?-A?) (12) 

u>S = \{A* i+ Af) (13) 

of the metric, the dreibein and the spin connection are made, and the metric is non-degenerate. 
Solutions with metrics that are degenerate somewhere need special study. In the present paper we 
encounter cases where the degeneration corresponds to coordinate singularity or to a conical singu- 
larity. In some of the cases the degeneration can be directly associated to horizons, with coordinate 
singularities in the accompanying 'Schwarzschild-likc' coordinate systems. Such degenerations may 
be handled by attaching another coordinate patch with a boundary and gluing them together by 
the appropriate matching conditions. Then one may find a new coordinate system covering the 
boundary region, with a metric which is non-degenerate. Thus the degeneration is not a coordinate 
invariant concept (unless restrictions are imposed on the allowed coordinate transformations at a 
supposed boundary of spacetime). 



2.2 Chern-Simons representation of the BTZ black hole 

Now we want to write down the Chern-Simon fields corresponding to the metric in each region, and 
then verify that the field strength F vanishes even at the horizons. We need F to vanish everywhere 
in the interior of our space except at singularities for the solutions to represent a spacetime with 
constant negative curvature. A non-vanishing field strength at the horizons can only come from a 
discontinuity in the A field when we glue the different regions together (recall that F = dA + A/\ A, 
and if A contains a step function the differential gives rise to a delta function). Since derivatives 
transverse to the boundary only appear in F for the longitudinal components A t and A^, it is 
enough to ensure that these components are continuous. 

Knowing the metric in the different regions, we may choose corresponding dreibeins and derive 
the corresponding spin connections from the equation of motions, de a + Lo a \, A e b = and dui a + 
(l/2)e a bc^ b /\lo c = —7rpe a bc eb A e c . The result is unique up to local Lorentz transformations, and 
a simple choice is 

j ( e — — sinh(/9 — a — ^)[r + dt — r_d</>]J + cosh(,o — a — §)[— T-dt + r + d<p]3x + rfpj 2 
y lo = — sinh(/9 — a — §)[— T-dt + r + d(j)]3o + cosh(p — a — ^)[r + dt — r_d</>]Ji 



II 



{e = df 
uj = si 



= dpJ - 
sin(p 



sin(p — a)[r_dt — 
- a)[r + dt — r_d(j>]Ji - 



i<(>]3i + cos(p — a)[r + dt — r_c?0]J 2 
cos(p — a)[r_dt — r + d0]J 2 



(14) 



III 



e = sinh(/? — a)[— T-dt + r + d(f)]J a + dpJi + cosh(/? — a)[r + dt — r_G?0]J 2 
lo = sinh(/9 — a)[r + dt — r-d<p]3 + cosh(p — a)[— T-dt + r + d(f)}J 2 



4 



These dreibeins and spin connections can be compared with those of Cangemi et al 27 , who use 
a different radial coordinate (the same as in the metric ([!])). Otherwise the differences are the 
choices of some of the signs and in the outer and inner regions the Lie algebra components are 
interchanged. Our choice of a means that closed timelike curves are excluded in the region p > 0, 
and it corresponds to the boundary condition (at p — 0) that the ^-component of the dreibein is 
lightlikc. In effect it relates the tangential components of A and A at these boundaries. 
From A = lj + e and A = to — e we get the Chern-Simons fields 

7T 7T 

I : A^ = A t = (r+ - r_)[cosh(/? - a - -)J 1 - sinh(p - a - ^)Jo] 

77 : A = A t = [r+ - r_)[cos(p - a)J 2 + sin(p - a)J x ] ( 15 ) 
III : Acf, = A t = (r + — r_)[cosh(p — a)3 2 + sinh(p — a) Jo] , 

and 

- - 7T 7T 

I : A ( j > = -A t = (r+ + r_)[cosh(/9 - a - -)J X + sinh(p - a - — ) J ] 

U:A^ = -A t = (r+ + r_)[cos(p - a)J 2 - sin(p - a)J x ] ( 16 ) 
III : A^ = —A t = (r + + r )[cosh(p — a)3 2 — sinh(p — a) Jo] • 

Here we see that our choice of dreibeins make the longitudinal components of A and A continuous 
when passing between the regions (/ — * II and so on). The A p and the A p just become 

I:A p = 3 2 II : A p = J III : A p = Jj 
I:A p = -J 2 II :A p = -Jo III:A P = -3 1 . 

Thus the only discontinuous component is A Pl which in fact can not contribute to the field strength 
since it only depends on p and the other components are continuous. F pp vanishes by antisymmetry 
and the off-diagonal terms F^ p and Ft p vanish by relating the discontinuities of d p A$ and d p At 
respectively with that of A p . 

There is an important distinction between how the boundaries between the regions are treated 
in the Chern-Simons formulation and in the metric formulation. In the metric formulation we 
can be forced to match regions with different forms of the metric (or find a coordinate patch 
covering the boundary) in order for Einstein's equations to make sense everywhere. A naive 
analytic continuation of the outer metric (||) to all p would divide spacetime in two disjoint pieces. 
In contrast, the Chern-Simons formulations seems to leave us with a choice. There is nothing 
wrong with the expressions for the vector potentials I, II or III, even if they are extended to all p. 
We can take those expressions as they are (giving us a problem in the gravitational interpretation) 
or we can match solutions and get the BTZ solution. 

From a Chern-Simons perspective the matched discontinuous solutions and the smooth solutions 
are indistuinguishable in the outer region, and they both make equally good sense in the interior. 
Only imposing boundary conditions in the interior or imposing special gauge conditions may pick 
out one solution as preferable to the other. Thus a sound gravitational interpretation of the 
solutions is only possible given special boundary conditions or gauge fixings of the vector potential. 
In generalizing the BTZ solution we will ensure that the boundaries of different regions are always 
matched in the same way as in this original BTZ solution. 

To prepare for more general solutions let us write the BTZ solution in cartesian coordinates, 
p = \J x 1 + y 2 and <f> = arctan (^). In the inner region we can write it as, 

A y =g^2+ 9 ^Ji (19) 
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where 



q = r+ - r_ (20) 

and 

g = [cosh( ( o — a)3-2 + sinh(p — a) Jo]- (21) 

The second vector potential 

A x = - d xP Ji (22) 

x z + y^ 

A y = g-^-d y pJ 1 (23) 
x l + y z 



where 



and 



q = r+ + r_ (24) 



g = [cosh(p — a)^2 — sinh(p — a) Jo]. (25) 

In cartesian coordinates it looks as if p = denotes a single point in space. There is no a priori 
justification for this since we chose p = to be special by hand, and all other equations p = const 
denote topological circles. On the other hand, we excluded p < on physical grounds, to get rid of 
closed timelike curves. Furthermore, calculating F in cartesian coordinates we get a delta function 
at the origin which we may formally regard as a source, and in this context we can also regard 
p = as a single point. 



2.3 Holonomies 

In a gauge theory of flat connections, F = F — 0, gauge invariant observables are scarce. The fields 
are locally pure gauge A = U~ 1 dU , for U an element of 50(1, 2), and any non-trivial observable 
has to be associated with the boundaries of spacetime or be topological in nature. The simplest 
topological observables are holonomies (or Wilson loops) measuring the effect of parallel transport 
along a closed loop in spacetime. For flat connections the result can only be non-zero if the loop 
C x (based at x) is non-contractible. Then the Wilson loop 

W(C X ) =Vexp(^(j> A^j =U- 1 (x)U(x + C x ) , (26) 

where V denotes path ordering of the exponential. As observed by Cangemi et al. [^Tj it is simplest 
in our case to take the closed curve C x at constant radial coordinate, i.e. along a level curve of p. 
For two curves C x and C v which can be continuously deformed into each other, but are based at two 
different points x and y, the holonomies are conjugate, W(C y ) = U(y)~ 1 U (x)W (C X )U (x)~ 1 U (y) . 
The eigenvalues of W for two curves which can be continuously deformed into each other are thus 
equal. These eigenvalues are determined by the parameters q, q and the eigenvalues of the SO(l, 2) 
Lie algebra elements g and g. It does not matter in which coordinate patch we follow the level 
curves, because we have ensured that the connections are flat also at the boundaries between the 
patches. 

Since the gauge group is a product of two rank one groups it is enough to characterize the 
eigenvalues by the two traces Tr W(C) and Tr W(C). For the BTZ solutions we obtain 

Tr W{C) = 2 cosh[7nj], Tr W{C) = 2 cosli^g], (27) 
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for Wilson loops in the two-dimensional representation of SO(l, 2). Via Equation || the holonomies 
are then related to the mass and spin of the black hole. In the complete classification of conju- 
gacy classes of SO(2, 2) Lie algebra elements ^ one finds that holonomies corresponding formally 
to imaginary q or q may occur, and furthermore that the case of coinciding eigenvalues (when 
Tr W(C) = 2 or Tr W(C) — 2) allows for non-trivial holonomy matrices (in addition to W = 1 
or W = 1). These cases can be dealt with in the Chern-Simons formulation by modifying the 
expressions for g and g. 



3 Multi-black hole solutions 

We will generalize the solution ( |l^ ) to the case were we have arbitrary many singularites. We will 
use the same form of the solution as in the inner region, regarding the Lie algebra direction of A.^ 
We may then try a solution 

A = dhli + (/ + dt)g (28) 
g = go(h)3 +g 2 {h)3 2 , (29) 

where h is a scalar function generalizing the radial coordinate p and / is a spatial one-form which 
is closed except at isolated sources 



N 
i=l 



df = 2ir}qiS 2 (x-Xi) dxAdy . (30) 



The qi determine the strength of the sources (the masses and spins of black holes). By integrating 
( |30| ) over a large disk D enclosing all sources we obtain 

N 

f= / d/ = 2fV 9 , = 2iQ . (31) 
dD Jd l=1 

If appropriate boundary conditions on / are assumed, / — > Qdcj) as p — > oo. Then we may regard 
A t = Qs(h) — ► A^ as a natural generalization of the relation A t — A^ — satisfied by single BTZ 
black holes. This is consistent with the ansatz ( ^8f ) after rescaling t. 

The equations of motion dA + A f\ A = arc satisfied by the vector potential ( ^8| ) outside the 
sources, x ^ Xj), provided 

-dh= 92 > -dh =9 °- (32) 

We recognize the equation for the hyperbolic functions entering the BTZ solution, but now their 
arguments have been generalized from p to h. By permuting the Lie algebra elements Ji in 
Equations (^8[ ^9|) one obtains solutions generalizing the BTZ solutions for all three regions, 
provided the signs in Equations (|2^, |32|) are changed accordingly. Matching of the regions 
works precisely as in the BTZ case. Note that the Lie algebra element g is spacelike, null or 
timelike depending on the sign of Tr g 2 , and that its sign is necessarily constant all over spacetimc 
for the present solutions. The one- form f/Q generalises the angular one- form dip in the BTZ case. 

Although any choices of h and of / satisfying Equation ( |30| ) are consistent with the equations 
of motion, we will concentrate on boundary conditions and combinations of A and A solutions that 
reduce to ordinary BTZ solutions both for asymptotically large h and close to the sources (regions 

3 If we exchange Ji with J2 we also have to change signs in front of the Jo component, in order to preserve the 
commutation relations which govern the equations of motion. If we exchange Ji with Jo we have to change the 
last condition in ( p2[ ) below. There will then be a minus sig n in front of gi, in effect exchanging trigonometric and 
hyperbolic functions. This is precisely the case in (^), ([lq) and (p^[). 
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of closed timelike curves). All the important new features of these generalized solutions are then 
associated with the fact that they are multi-centered, which in its turn implies that there will be 
critical points of h and /. Such critical points can give rise to degenerate metrics, a subject we 



shall return to in Section (4.2) 



The second gauge field A has analogous solutions in terms of h, 52 (ft) and go(h). In order to get 
solutions similar to the BTZ solutions we may choose h = —h, 52 (ft) = 52 (ft) and 50 (ft) — <?o(ft)i 
guided by Equations ([15]) and (|l6|). In an inner region with g 2 and go even and odd functions 
respectively, we obtain the vector potentials 



A = {f + Qdt)g {h)3 + dhJ 1 + (/ + Qdt)g 2 {h)3 2 
A = -{f- Qdt)go(h)J - dftJi + (f - Qdt)g 2 (h)3 2 



(33) 



and the metric, 

ds 2 = -g 2 (h(x, y)) {r^dt - f+(x, y)} 2 + g 2 2 (h(x, y)) {r+dt - f^(x, y)} 2 + dh(x, y) 2 , (34) 
where 

f±f Q±Q 
f± = —> r± = —2~ ( 35 ) 

The metric is easily compared with the BTZ metric (^|) in the inner region (III). The function 
h + a has replaced the radial coordinate p, go and g 2 represent the hyperbolic functions, and r±d<p 
is replaced by f±. The last change is the most significant one, since two different one- forms are 
needed to generalize dip. Only when /+ and /_ are proportional do we get a direct multi-source 
generalization of dep. This happens when the ratio of the two charges at each source is constant. 
Irrespective of this we can make direct contact with the BTZ-solution very close to a charge, where 
the effect of the other charges is negligable, or at asymptotically large distances, where the sum of 
the charges dominate the solution. 

In the general case ( |34j ) we can still define regions of type /, II and III, between which the 
solutions have to be matched, and different choices of the function h gives different regions (even 
their topologies may be different), but they are actually related by gauge transformations, as we 
proceed to discuss. 

3.1 Gauge transformations 

One can check that the gauge transformation 

5 A = dShJx + [A, ShJi] (36) 

amounts to a change of h into h + 8h in the solution for region implying that solutions with 
different functions h are equivalent if only their boundary conditions are the same. Of course, 
an analogous statement is true for A. We stress that the solutions are only equivalent in the 
Chern-Simons formulation of pure gravity. To see this we may study horizons. 

The boundary between region / and II resembles a horizon and it is actually an event horizon 
for constant charge ratios of all sources, if it consists of a single connected component. This is 
because we can then use coordinates h and ip with dtp = f to obtain the ordinary BTZ metric in 
the exterior region. The transformation to these coordinates works asymptotically and also in the 
whole exterior region provided / does not have a zero there. In fact, we will show later that the 
multi-black hole solutions have singularities at zeroes of /. These singularities may be inside or 
outside a physical event horizon depending on the choice of the function h. Such a difference could 
for instance be detected by the propagation of light rays in the background metric. Light rays are 
of course not included in a Chern-Simons description. 

Even if there is little physics in the function h, the multi-black hole solution also depends on the 
forms /+ and /_, which in their turn depend on the positions and charges of the sources. As will be 
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discussed in the next subsection, the charges may be directly measured by holonomies around the 
sources. The positions of the sources are trickier, and cannot be resolved by the holonomies. Other 
available observables are the asymptotic charges |^8| . The general /+ and /_ are asymptotic to the 
corresponding BTZ forms, and the issue is if the approach is fast enough to give finite asymptotic 
charges, but also slow enough to give non-zero values. 

As an example we may compare a single source BTZ solution A\ with a solution A 2 with 
sources separated by a small coordinate distance xo in the x direction. Then one finds 

A 2 = A 1 + S^Ai =A X + x (q 2 - qi)d [ V ) g . (37) 

\x £ + y J 

Thus <5i2j4i0 scales as gp -2 while Ai^ scales as g/f -1 with p and the change is subleading. If the 
change S^Ai^ can be written as an infinitesimal gauge transform 6\ l2 Ai with a decreasing gauge 
parameter A then the separation of the two sources is truly a matter of gauge choice at infinity 
and it is not detectable by any asymptotic charges. (It will still be detectable by holonomies, 
corresponding to the fact that the gauge transformations are not defined everywhere, or do not 
belong to the identity component of the gauge group.) The problem in our case is that the 
asymptotic behaviour of the BTZ solution implies that g has an exponential dependence on p. The 
same is true for A12. Then the boundary values of the fields and the transformation parameters 
are not well defined. Fortunately, this problem may be circumvented by discussing the vector 
potentials 

A' = e pj2 de- pj2 + e pj2 Ae- pj2 , (38) 

which locally are gauge transforms of A but satisfy different boundary conditions. In fact A' BTZ 
is a constant and the A' of our generalized multi-source solutions approach constants at infinity. 
The A' do however give rise to metrics which are everywhere degenerate, and we just regard 
them as auxiliary solutions which help distinguishing asymptotic gauge transformations and global 
transformations generated by asymptotic charges. The parameters of global transformations on A' 
go to constants at infinity while true gauge transformations vanish asymptotically. The effect of 
both kinds of transformations on the fields A is simply obtained by the mapping inverse to ( p8| ) . 
Conversely, by mapping to A' transformations on A may be classified as gauge transformations or 
global transformations (or as changing boundary conditions). 

Returning to <5i2^4i0, its image S^A'^ under the map ([58]) vanishes at infinity, implying that 
asymptotic charges are left invariant by moving sources apart. In fact, S^A'^ = S\ l2 Ai for 
A12 = —yg'/ {x 2 +y 2 ) with a constant g'. Since A i2 diverges at the origin it does not give a globally 
well defined infinitesimal gauge transformation and there can still be a physical difference between 
the solutions. In conclusion, solutions with different numbers of sources are inequivalent because 
of different holonomies, while different positions of the sources may or may not be observable 
depending on the global properties and boundary conditions of the finite gauge transformations 
effecting the translations. The asymptotic charges are insensitive to these details, so they may 
be thought of as generating transformations common to several different sectors labeled by the 
numbers of sources, and possibly by their positions. 

3.2 Multi-black hole holonomies 

We now wish to calculate holonomies 



c 



Tr V exp / A (39) 



We first calculate the ordinary integral over a closed loop, 

d h(x,y) \ ( dh{x, y) 
—3 1 jdx+\f v8 +— Bjj 



A= I [(f x S+^^^)dx+(f yg +^^3 1 )dy) , (40) 
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here written out for an inner-type region. If the function h(x,y) is choosen in such a way that 
there are closed level curves of h(x,y) the term dh in the integral is zero, and furthermore g is 
constant. Then the integral depends on which charges qi are enclosed by the level curve, 

/ A = g f f = 27Tg y; Qi = 2n S q c , (41) 
Jc Jc ig/c 

where Ic denotes the set of enclosed sources and qc the enclosed charge. Since the eigenvalues of 
the traceless real matrix g are necessarily both real or both imaginary and add up to zero, we may 
write 

Tr (Vexp J A^j = e 27rqcX + e~ 27rqcX = 2 cosh (27rq c X) (42) 

where A is one of the eigenvalues of g, and independent of h. The matrix corresponding to the 
A must also have either both imaginary or both real eigenvalues which we call A and —A. So in 
general we get three different holonomy types depending on the eigenvalues A and A: either one is 
real and one imaginary, both are real or both are imaginary. When we just have one singularity it 
is known that these types will correspond to different quotients of anti-de Sitter space. Bahados et 
al have shown how different spaces are obtained from anti-de Sitter by modding out subgroups 
of SO(2, 2), and that BTZ black holes belong to one of these classes of spaces. They also find three 
different types of spaces. The correspondence between their eigenvalues A' and our eigenvalues is 
A^ = qX — qX and X' 2 = qX + qX. In our language the generic BTZ black hole corresponds to the 
case with two real eigenvalues. When both are imaginary we generally get conical singularities, 
except in the case qX = qX — i/2, which curiously corresponds to AdS space^j. In fact, we may also 
find 'multi-AdS solutions' with several of these AdS charges. They may possibly serve as ground 
states of multi-black hole sectors. Note that the holonomy around a single AdS charge is almost 
trivial, and around two it is entirely trivial. 

Notice that we have not mixed holonomy type for the different singularities. It would be 
interesting to find solutions where the sources give rise to different types of holonomies. 



4 Two sources 



We will study the solutions for the case with two sources in more detail. After verifying that 
the solutions approach the single-source solution asymptotically and for vanishing separation of 
the charges, we will continue with a generalization to several sources of the procedure to exclude 
CTCs, and we will also discuss how the multi-source solutions generically contain additional (mild) 
singularities. 

Solutions with sources at x = x\ — xq and x — X2 — —xq can be written, 



At = [(r 1+ - n_) + (r 2+ - r 2 _)]g (43) 
A x = (h x + f 2x ) S + ^^J 1 (44) 

A y = (f ly + h y )g+^^J 1 (45) 



where 



f ~y f x-x 

= qi i \2~r^i hv = 91 ? \2 i 2 ( 46 ) 

(X — Xo) +y [X — Xq Y + 2T 



J 2* - & , „^ 2 , „,2 fa - 92, \ 2 - 2 ■ (47) 



-_y , _ x + xp 

' (x + xq) 2 + y 2 y (x + x ) 2 + y 

4 M = — 1 and J = are obtained from Equations |24| ) and (^|) and the metric (Q) then represents AdS 
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The conjugate field A, 

^ = (/i* + /2*)g + ^jg^Ji (48) 
A y = (h v + hy)g+^^-Ji (49) 

where, 

g = 92{h{x, y))3 2 + 9o{h{x, y))J . (50) 

In the BTZ-like inner region with h = —h, g^Qi) = g 2 (h), go(h) = go(h) and g 2 and go even and 
odd functions respectively, we find the metric 

ds2 = - go( h (x,y)) {(r x _ + r 2 -)dt - {r 1+ f lx + r 2+ f2x)dx - (r 1+ f ly + r 2 +f2 V )dy} 2 
+ g 2 ( h ( x ,y)) + r 2+ )dt - {r X -f lx +r 2 _f 2x )dx - {rx-fxy + r^f^dy} 2 

dh(x,y) dh(x,y) ^ ' 

— a dx H a dy 

ox ox 

So far the function h has been left unspecified. If, for instance, we choose h(x,y) — y 'pip 2 — a 
in terms of the radial coordinates p\ and p 2 centered on each of the two sources and a function 
a approaching a constant (^) at infinity and in the limit xq — > 0, we can ensure that the BTZ 
solution is approached both at infinity and as xq — ► 0. To verify this, start by looking at the metric 
in the outer region 

ds 2 = - sinh 2 ( v / / dip 2 - a) ((n+ + r 2+ )dt - ({rx-hx + r 2 -f 2x )dx + (ri-fx y + r 2 -f 2y )dy)) 2 

+ cosh 2 (^/p 1 p 2 - a) ((n_ + r 2 -)dt - ((r 1+ f lx + r 2+ f 2x )dx + (n+fiy + r 2+ f 2y )dy)) 2 ^ 

( dh(x,y) dh(x,y) 
+ {^x- dx+ ^x- dy 

to see how it behaves asymptotically at infinity. In terms of polar coordinates (p, </>) centred around 
(x,y) = (xi,0), (implying p = pi) 



_ —y _ —psmq) _ x — x _ p cos <p 

fix — o" — 5 Jly — 5 — 2 

pZ pZ pZ pZ 

—y —psiiKp X + Xq pCOS(j) + 2xo 

J2x = 0" — 2 J 2 V = 2 = 2 

p 2 p 2 p 2 pi 

the metric takes the form 

2 2 , , \\ fr 2 -(p 2 +2x p cos <p) +ri-pl 2r 2 -X sm^ 

ds outer = — smn Vv P/°2 — oc) < r+dt — I ^ dq> H j dp 

I V P2 P2 

+ cosl.\^ 2 -a)^dt- ^ 2+{p2 + 2X0P p^ +ri+pl 



(53) 



f dh&p) | dh{cj>,p) 
\ dcj> dp 



dp 



(54) 



We see that the metric is asymptotic to the BTZ solution with r+ = ri+ +r 2 + and r_ = n_ +r 2 - 
when p — > 00 or so — > 0. 
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4.1 Exclusion of closed timelike curves 



In the BTZ solution (||) there are closed timelike curves for p < 0, and we expect similar pathologies 
in the multi-black hole solutions inside the black holes. It is natural to cut off the range of the 
coordinates precisely where CTCs are encountered. Here we show how this can be done in the 
case of two sources. The same procedure can be used for any number of sources. The resulting 
spacetimes then have singularities in the causal structure if they are continued 'inside' the sources. 

Just as for the BTZ case (||) we need the vector field d<j, for some periodic coordinate <fi to become 
lightlike at each source in order to exclude regions containing closed timclikc curves. Coordinates 
which are periodic around curves enclosing only single sources are readily found. We may use the 
angle between the line from the source to a point and the positive x direction, or we may use df+ 
and df- to measure angular differences. Close to the sources these measures of angle all agree up 
to proportionality constants. 

To localize the causal singularities to the positions of the sources it is then enough to choose 
the function a appropriately. In order to encounter closed timclikc curves we have to go to the 
inner region. 

First study the metric in the inner region. It is obtained from the outer metric (Q) by ex- 
changing r + with r_ : 



ds 2 

""inner 



sinh 2 ( A /pp 2 — a) < r_di — 



+ cosh 2 ( v / pp2 — a) yr + dt ~ 
+ {-^ d ^ + ^p— dp 



r2+{p 2 + 2x p cos 4>) +r 1+ pl 



r2-{p 2 + 2x pcos(/)) + ri-pl 



pi 



2r 2 +x sin <p 

2t2-Xq sin tfi 




(55) 



Now take a look at the component, 



r2+ (p 2 + 2x p cos 4>) + rx + p% 



Pi 



Pi 



7*2- (p 2 + 2x a p cos 4>) + ri_/j 2 



sinh 2 (y/pp 2 — a) 



cosh 2 {\fp~p~2 — a) 



(56) 



£Qpsm 2 20 /<9cA ^x%^fpsva.2(f> f da 



Pi 



3/2 
P2 



when p = 0. We must choose a in order to make the vector field lightlikc at x = xq. For a with 
a = when p = 0, the condition that d^, becomes lightlikc becomes 



a = arctanh 



(ri + ) 2 sinh 2 a + (ri_) 2 cosh 2 a = 



(57) 



In the same way we can change to polar coordinates centred around x — —Xq which instead would 
lead us to the condition, 



9U 



'{ r 2+) 2 sinh 2 a + {r 2 -) 2 cosh 2 a = 



a = arctanh ( - 2 — 
V r 2+ 



(58) 
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Figure 1: The inner and outer 'horizons' in the xy-plane at fixed t for different xq. 



In order to have both these conditions satisfied a can only be a constant in the case r\—/r\+ = 
7"2- l r 2+- Still, there are many ways of choosing an a(p, <fi) that does not affect the singularities or 
the asymptotics of the solutions. We may choose a to be a constant at infinity, for instance 

a = arc tanh ( ri ~<* ± ^ ) (59) 

We see that in the case r\-/r\+ = r2-/r 2 + this a will reduce to a constant. This will also be the 
case when xq = 0, i.e. when the singularities are in the same point. The requirement d^, a = 
when p = is also easily seen to be fulfilled. 

To make the analogy with the BTZ case complete the different regions we had can be generalized 

to, 

I : < p < a < sjp\pi < a 

II : a < p < a + § => a < y/pTpl. < a + f (60) 
III : a + | < p ^a + f< ^/pTpl 

In figure [l] we have plotted the 'horizons' when we have fixed r+ and r_ but varying distances xq 
between the singularities. Although the equations determining the boundaries of the regions are 
similar to the single-BTZ case we cannot be certain that we are dealing with true horizons, unless 
we trace light rays through the new geometries. This explains the quotation marks. 

4.2 Singularities 

The metric ( |34| ) may locally be written 

ds 2 = -gldT 2 + gldti 1 + dh 2 , (61) 

with dT = r-dt — /+ and d$ = r + dt — /_, since / + and /_ are closed forms. This metric 
degenerates where go or g 2 vanishes, where one of the functions T(t,x,y), &(t,x,y) or h(t,x,y) 
has critical point, and where dT, d$ and dh are linearly dependent. The coordinate singularities 
at the BTZ horizons and their multi-black hole generalizations belong to the first case, but our 
solutions also display the other types of degeneracies, and we now proceed to investigate their 
interpretation. 

In the case when one of the functions T, $ or h has a critical point, one may ignore the effects 
of the functions go or g 2 locally, since they may be absorbed into redefinitions of T, $ or h only in 
exceptional cases at the expense of changing the nature of the critical point (but see the following 
subsection to appreciate the importance of these exceptions!). Then the singularity is precisely of 
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the kind discussed by Horowitz |£9[ for zero cosmological constant. The simplest such singularity 
occurs between two equal charges separated by some distance. 

To see what happens we study the equal charge solutions close to the origin. There /+ = 
/_ = because the contributions from the two charges cancel by symmetry. The metric ( |5lj ) 
then degenerates at the origin at all times, because dT and d$ both become parallel to the Killing 
direction dt. Furthermore, h, which approaches infinity at infinity and assumes local minima at the 
positions of the charges, has to have a saddle point. Due to gauge invariance (^) the position of the 
saddle point may be chosen to be at the origin, making the metric on this line (in spacetime) even 
more degenerate, of rank one. Generically we instead expect degenerations to rank-two metrics 
on two-dimensional surfaces [^9|, [30|. In fact, we have found that the map (t, x, y) — ► (T, h) has 
three singular fold surfaces joined pairwise at three cusp lines if the saddle point of h is displaced 
slightly. The geometries of such complicated singularities deserve a special study, but for our 
purposes it is enough to find the simplest singularities in a gauge equivalence class. 

Returning to the case of coinciding saddles we proceed to determine the geometry close to the 
saddles. There we have approximately 

h =ax 2 — by 2 

, v (62) 
r_/ + =r + /_ =cd{xy) . 

By rescaling coordinates and h we find a spatial line clement 

ds 2 = d(xy) 2 + - A d{x 2 - y 2 ) 2 = (x 2 + y 2 ) (dx 2 + dy 2 ) (63) 

The area Ao and circumference Co of circles around the origin are then related by Cq = 8irAo in 
contrast to the Euclidean relation C 2 = AttA. Since the metric is manifestly flat the difference can 
only be due to a conical singularity at the origin, and we conclude that there is a negative deficit 
angle of 2tt. 

We have argued that simple conical singularities with a surplus angle of 2n appear in the 
geometries with two equal sources provided the gauge is chosen so that saddles of h coincide with 
zeroes of /+ and /_ . For n sources h typically has n — 1 saddles since it is chosen to have n local 
minima at the sources and a maximum (infinity) at infinity. Similarly / + and /_ typically have n— 1 
zeroes, because of the n sources and the behaviour at infinity. If /+ and /_ are proportional their 
zeroes coincide, and h may be chosen to have saddles at the same points. Fixing the behaviour of h 
appropriately close to its saddles the local calculation is then the same as between two sources, and 
we conclude that there are n — 1 conical singularities. Physically the proportionality of /+ and /_ 
means that the sources all have the same ratio J/M of spin and mass. Other source distributions 
generally lead to more complicated singularities in the geometry. Some of these may be removable 
like the coordinate singularites of the BTZ geometry, but some are likely to be required by global 
arguments, like the conical singularities we have just discussed. 



4.2.1 Geodesic singularities 

The stationary conical singularities discussed above have been found before by Clement and 
by Coussaert and Henneaux p(| . These authors have also remarked that such singularities do not 
follow geodesies. This is quite disturbing for the commonly used hypothesis that Chern-Simons 
theory should be relevant to the counting of black hole states. Already at the classical level would 
Chern-Simons theory give rise to geometries which seem to leak energy and momentum! 

Fortunately, Chern-Simons theory itself contains the answer to the problem. By asking under 
what precise conditions the singularities are non-geodesic we may find an exception: when the 
singularity is located at a horizon. This case was already mentioned by Clement, but not in 
the Chern-Simons context where it becomes truly important. While the set of geometries with 
singularities fixed to horizons may seem like an exceptional set of measure zero, in Chern-Simons 
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theory they are not exceptional. In fact, large class of solutions (and all those considered by 
Coussaert and Henneaux) may be written in a gauge such that the singularities are located at the 
horizon and thus follow geodesies. We now proceed to give some details of this argument. 

We need to evaluate the Christoffel symbols Tf t and Y v tt which vanish precisely where there 
are static geodesies. If they can be made to vanish at the conical singularities the puzzle of the 
unphysical geometries is solved. We study the Coussaert-Henneaux solutions, which are essen- 
tially ordinary BTZ solutions, but with the mass and angular momentum distributed in the same 
proportions on several sources. In our language this means that 

r-.f+ = r+f- = r-r+f , (64) 

where the single form / encodes the source distribution. As has been pointed out several times 
above this assumption simplifies the interpretation of the solutions considerably. Now 

T x - ±(r 2 - r 2 ) (65) 

where the sign (and the label i) depends on the region. In general this expression and the one 
for Ijj diverge at a common zero of / and critical point of h, but if g or g' vanishes at the same 
point the whole expression instead goes to zero. This is what happens if the conical singularity is 
located at a horizon. 

It only remains to argue that the singularities can be moved to a horizon. Indeed, in region III 
a infinitesimal shift of h is equivalent to an infinitesimal gauge transformation (^), and similar 
relations exist in the other regions. Assuming that these transformations can be integrated, we 
conclude that changes of function h are gauge transformations. By adjusting h we can then make 
g or g' vanish at a conical singularity, i.e. a gauge transformation may take the singularity to a 
horizon, where it follows a (null) geodesic simply by being stationary. 



5 Conclusions 

We have constructed and investigated solutions to three-dimensional AdS gravity which generalize 
the BTZ solution. While the ordinary BTZ black hole can be viewed as a single source solution in 
the Chern-Simons formulation, we have constructed multi-source solutions. These solutions give 
rise to a kind of multi-black hole solutions, which however also display other singularities. In the 
simplest cases the additional singularities are fixed conical singularities, but more complicated cases 
also occur. Einstein's equations break down at these singularities, so they represent geometries 
which are not allowed in pure einstcinian gravity. On the other hand, they occur very naturally in 
the Chern-Simons framework, which is natural for quantization, so we believe that these multi-black 
hole solutions should be included in a full Chern-Simons treatment of BTZ black hole entropy. 

We have also shown that a large class of these multi-black hole solutions allow a gauge choice 
which ensures that the singularities in the corresponding geometries follow geodesies. Geometrically 
the solutions then precisely encode the BTZ solution outside a number of horizons. These horizons 
are however all connected with each other, since the conical singularities which join them can 
not appear outside the horizons without violating the geodesic equation. The union of all these 
horizons appears to the outside observer as a single horizon. Only at the horizon (and inside) 
is the difference to the single black hole solution noticable. In this picture of a single horizon, 
special light-like geodesies on the horizon are identified pairwise, since they in fact represent the 
same conical singularity, only approached from two different directions (two different ridges on the 
saddle point of the function h). 

Although we have not attempted in this paper to find the quantum states corresponding to the 
multi-black hole solutions, we have provided evidence that such states should be included in the 
black hole spectrum. Namely, the asymptotics at infinity of the classical solutions approach the 
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single-BTZ solutions so rapidly that the difference can not be detected by any asymptotic charges. 
Only non-asymptotic observables like the holonomies distinguish between the solutions. It then 
seems quite unnatural to exclude the sectors with multiple sources, in particular since the sources 
may be hidden inside the horizon. Presumably, the additional sectors of the boundary conformal 
field theory that are required to represent multi-black hole solutions can also be understood by 
purely two-dimensional considerations, for instance by the requirement of modular invariance. 
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